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[ We consider functionals of long-range dependent Gaussian se- 

>0 . quences with infinite variance and obtain nonstandard limit theo- 

rems. When the long-range dependence is strong enough, the limit 
is a Hermite process, while for weaker long-range dependence, the 
limit is a-stable Levy motion. For the critical value of the long-range 
dependence parameter, the limit is a sum of a Hermite process and 
a-stable Levy motion. 









1. Introduction. Define the stationary Gaussian sequence 



> ; 

OO I wliere the ^i are independent standard Gaussians and bj = j^~^''^Li{j), 

J^ ' with ^ < H < 1, Li{i) slowly varying and X^j^o^j — ^- ^^^ {-^i) ^'^^ then 

«vj ■ long-range dependent with Hurst index H. Denote the filtration generated 

■• I by {^i)i£Z as (J-'i). The Hermite polynomials form an orthogonal basis for 

Q ■ L^(M,e~^ '^), so when / is a function such that Ef{Xi)'^ < oo, the chaos 

OO , decomposition of f{Xi) is given by 
O 



(1.1) f{Xi) = Y,fkhu{Xi), 

k=0 



j^ ■ where the sum is interpreted by convergence in L^ and fk = ^E{f{Xi)hk{Xi)). 

The Hermite rank of / is the smallest k>l such that /fc 7^ 0. For more details 
on Hermite polynomials and their relationship to Gaussian Hilbert spaces, 
see [8]. The limits of normalized sums of the form J2i=i fi-^i) ~ ^fi^i) were 
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established in [5, 14] and [2], when Ef{Xi)'^ < co, in terms of the Hermite 
rank k and H . When 1 — k;(1 — H) > 1/2, the hmit is a Hermite process, but 
when 1 — «;(1 — H) < 1/2, the hmit is standard Brownian motion. 

This paper considers the limits of normalized sums when the f{Xi) have 
power-tailed distributions with index < a < 2. These functions / satisfy 

P(/(X)>x)~i±^L2(x)x-°, 
(1.2) _ 

P{f{X) < -x) ~ i^L2(x)x-" 

as x — > oo, L2 slowly varying and with /3 G [—1,1]. We can find constants a„ = 
n^'°'L^{n) such that P{\f{X)\ > an) ~ n~^ and L3 is slowly varying. We will 
focus on the case when 1 < a < 2 and / is centered so that E{f{X)) = 0. 

For background, a random variable X is stable if cX = Xi + X2 for some 
constant c and where Xi and X2 are independent copies of X. Such distri- 
butions are completely categorized and the non-Gaussian stable laws have 
characteristic function 



Eexp{itX) 




TTCa \ 

1 — «/3(sign(0)) tan — \+i^9 , a 7^1, 

-a°|0r(l + ip{sign{e)) \n\e\) + i^iO], a = 1, 

for 0<a<2, /?G [—1,1], a > and fiGM. We will denote this distribution 
by Sa{a,[5,iJi), following [11]. A process is an a-stahle Levy motion if it is a 
Levy process with increments which are stable; see [11] for more details. 

The standard approach which relies on the chaos decomposition does not 
apply here. We use the hypercontractivity of Gaussian Hilbert spaces to 
establish that the extreme values of such processes are asymptotically inde- 
pendent. The limit can be either a Hermite process, a-stable Levy motion 
or a sum of both. This striking result, where normed sums of power tailed 
variables with infinite variance may converge to a self-similar process with 
finite moments of all orders, highlights the complexity associated with the 
domains of attraction of the Hermite processes, all of which have this prop- 
erty. Little is known about the domains of attraction beyond the case of 
fractional Brownian motion (k = 1). Note that, in the case of i.i.d. variables 
(the classical central limit theorem), the only self-similar processes which 
appear as limits are Brownian motion and the a-stable Levy motions for 
0<a<2. 

Both infinite variance power tails and long-range dependence lead to par- 
tial sums increasing faster than 0{y/n). The different limit processes reflect 
the relative importance of these two effects. The composite process occurs 
only when the effects are balanced. 

When E\f[X)\^ < 00 for some p > 1, we can still define the fk and equa- 
tion (1.1) can be interpreted as a stochastic distribution; see [7], Corollary 
2.3.8. We give the Hermite rank its natural extension to L'^. 



oo, then 
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Theorem 1.1. Ifl-K{l-H)>^orifl-K{l-H) = ^andlim. ^^^^^ 



lnt\ 

1=1 

where Rk,h is the nth Hermite process given by the multiple stochastic inte- 
gral 

RK,H(t)= I I f\l{s>Xi)\s-Xif-^''^dsdB{xi)---dB{x^), 

JR« JO ~J^ 

B{t) being standard Brownian motion. If I — k(1 — H) < — , or if 1 — k(1 — 
H) = ^ and lim„ ^^^ = 0, then 

[ntj 



where R*{t) is a-stable Levy motion with R*{1) = Sa{{— — ^^^^^^^^^)^'",/3,0). 

in) 



Finally, if I — k(1 — H) = - and lim„ ^ \^\ = A G (0, oo), then 



[nt] 

a-1 Y. fi^i) ^ A/«i?,,H(t) + R*{t), 

where Rk,h o-nd R*{t) are independent. Convergence is taken to mean weak 
convergence on D[0,T] in the Skorohod Ji-topology. 

One might have expected the process f{Xi) to be in the domain of attrac- 
tion of a fractional stable process. However, it does not exhibit clustering of 
extreme values, as can be seen by a calculation of its extremal index. Our 
results imply that 



lim P 



max f{Xi)> a^n^j 



l<i<n' 



which shows that the extremal index is 1. On the other hand, a continuous 
version of Proposition 2.1 from [3] implies that fractional stable motion has 
extremal index strictly less than 1. 

The situation when < a < 1 is simpler as only convergence to a-stable 
Levy motion is possible. While not explicitly mentioned, the one-dimensional 
version of the case < a < 1 can be shown to follow from Lemma 5 and 6 
of [4]. 
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Theorem 1.2. IfO<a<l,then 

[nij 



«n'E/(^0^^*W' 



while if a = 1, then 

[nt] 



z-1 ^[/(X,) - E{f{X,)I{\f{X,)\ < a„))] - ^ - R*{t), 

TT 



= 1 



w/iere ■0 = In vr + /q°° u "^ {saiu — ul{u < 1)) du . Again, R*{t) is a -stable Levy 

motion with R*(l) = 5a(( ^^^"°i!°^^""^^^ )^/°,/3,0) w/ien < a < 1 and i?*(l) = 
5i(7r/2,/3,0) w/ien a = l and, again, convergence is in finite- dimensional 
distributions. 

The situation when < H < -^ is simpler with convergence to the stable 
limit in Theorems 1.1 and 1.2. 

The results herein remain true when we take any Gaussian sequence in 
the domain of attraction of fractional Brownian motion. However, the proofs 
involve much more tedious technical details and are no more informative, so 
we have confined our attention to the slightly less general result. 

There are other results in the literature where functionals of long-range 
dependent processes have been shown to have both Gaussian and a-stable 
Levy limits depending on the parameter of long-range dependence. For in- 
stance, Surgailis (in [13]) found this behavior for certain bounded functionals 
of long-range dependent moving averages of heavy-tailed random variables 
and (in [12]) for the empirical process of another moving average of heavy- 
tailed random variables. 

As a simple example, consider the case f{x) = \x\^' — £'(|Xj|'"). Proposi- 
tion 3 of [6] showed that ii j < H <1 and r > 0, then 

(1.3) ni-2^ Y. /(^O ^ f2R2,2H-l{l) 

i=\ 

as n ^> oo. With our results, we can extend this to all r. By Theorem 1.1, 
equation (1.3) holds when 1/(1 — 2i/) < r < 0, but 

-) n^-^^5:/(X,)^i2*(t) 



when -1 < r < 1/(1 - IH) and 

\nt\ 

Y^fiXi)^^-] R*{t) + f2R2,2H-l{l) 



n 

i=l 
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when r = 1/(1 — 2H). Also, using Theorem 1.2, when r = — 1, 



vr. 



lnt\ 



and when r < — 1, 



XA 



s(ix,rifix,r< (| 



r/2 



n 






9 ^ T-/2 L™*J 

-I ^^Ti^^^i 



vr 



CR*{t). 



2. Proof. Fix a d such that Y.f=db] = < min{a, 1/F} - 1. Then, by 
Theorem 5.1 of [8], the map f{Xi) i— > E{f{Xi)\J^i-d) is a hypercontraction 
from Li+^ to L^, so E(E(/(Xi)|^,_rf))2 < oo. With f/, = 0-1/2 ^~^6,-^,_,-, 
and expanding hk{Xi) according to Lemma D.l of [7], we have E{f{Xi)\J^i^d) 
EkLofkO'^^'^hkiUi). Then, by Theorem 5.6 of [14], when I-k(I-F) > 1/2, 



[ntj 



(2.1) L(n)-«n-(i-''(i"^)) ^ E{f{X,)\J^,_,) -. f^R.^t) 



i=l 



weakly as n — > oo. 

For our asymptotic independence result, we map our Gaussian sequence 
according to 



X.=Xr^=Y.ny%AB 



i + 1 



n 



B 



■J 



n 



where B{t) is some fixed Brownian motion. Brownian scaling guarantees 
that the distributions of the Xj's do not depend on n. With this definition, 
it follows from Lemma 4.5 of [14] that the convergence in equation (2.1) can 
be taken as convergence in L^ pointwise in t. We will routinely suppress the 
dependence on n of various objects. 

Fix a c > and let i/„ be the simple point process on (R\(— c,c)) x M 
given by point masses at points {f{Xi)/an,t/n). Let z/ be a Poisson point 
process with parameter measure 



i''{dx,dy) 



a^—\y\ dxdy, y > 0, 

1 — d 
a — - — \y\~°'~^ dxdy, y<0. 



Lemma 2.1. Let A be a finite union of finite intervals in M.\{—c,c) 
and let Zi = I[a~^f{Xi) 6 A]. Let V be a finite variance random variable 
measurable in the a-algebra generated by B{t). Then, 



(2.2) 



E 



E 



\nt\ 

Vj2{Zi-EZ^) 

i=l 



^0 
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as n —> oo and 

[nt] [ntj 

(2.3) EY^Y. ZiZ.^yif 

i=l jf=i+l 

as n^ oo, where fi/\ = i^'(A x [0, 1)). 

Proof. It is sufficient to prove the result when V G H \ the homoge- 
neous chaos of order / of the Gaussian Hilbert space generated by B{t) (see 
[8]). By equation (1.2), 

(2.4) nEZi -^ ^A 

as n — > oo. Zi can be written as J2k^o dkhki^i) ■ Applying [8], Theorem 
5.1 with the map Xi i-^ (lnn)~'^Xi, we then obtain X]iSo(l^'^)~^'^9fc^' — 
ll^i|l?+(inn)-2 = {EZi)^/(^+^'^'^r'^ = 0(n-2) and so glkl = 0{{lnn)^^n-^). 

Hence, since -B(E1=i ^fc(^i))^ is 0(L(n)2'^n2(i-'^(i-^))) when 1 - k(1 - 
H) > 1/2 and 0{n) when 1 - k{1 - H) < 1/2, we have 

/[ntj p \2 p Ant] \2 

E(YJ29khk{X^)) =Y9kE[Yhk{Xi)\ 

\i=l k=l ) fc=l \i=l / 

= C>(ln(n)2L(n)2'^n-2'^(^--^)) 

= 0(1). 

We choose p large enough so that (p — /)(1 — H) > 3. Now, decompose Xi 

as Xi = jiUi + 6iWi, where EU^ = EW^ = 1 and -^Ui = E{Xi\J'Q). So, 



E[vYgkhk{X,)T^=E[vYgkYi^^l^~'Hhk^jmhi{W,) 

\ k=p ) \ k=p j=0 ^-^ ^ 



^0 



Y.9kY.[A^^~'^ihk-m)E{Vh,{W.)\To). 
k=p j=o ^•' '' 



By [8], Theorem 5.10, EiVhjiWi))'^ is uniformly bounded over i and j. Since 



%' = E?^i 6? ~ (2 - 2H)f^-Hi{tf, we have 



sup I . ) 7j ■^ = o(n -^) 

and so when i > n^'^, 

/ oo ^ / T \ \ 



n^ 



fc=p j=0 ^ -^ ^ / 
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And, since EZf = 0{n~^), it follows that 



^0 







EiEivJ2gkhk{Xi) 

\ \ k=p / / 

as n — > oo, which proves equation (2.2). By [8], Corollary 5.7, 

and when \i - j\ is large, EX^Xj > 0, so E[ZiZj] = Y.Zo9k^KEXiXjf > 
Qq = [EZiY, from which equation (2.3) follows. D 



Proposition 2.1. Let A* he a finite union of finite rectangles in (]R\(— c, 
c)) X M. Let V he a random variable, measurable in the a-algebra generated 
by B{t) with \V\<1. Then, 

E[VI{vn{^*) = 0)] ^ EVE[I{v{^*) = 0)]. 

Proof. For e > 0, we can partition A* = IJ^^^ A^ = IJ^^^ Aj x [sj,tj), 
where Aj are finite unions of finite intervals in M\(— c, c), si < ti < S2 < 
■ ■ ■ <tk and J2j=i |/^A fe ~ ^j]"^ < ^- ^^ ^^^ write 



L{un{A*) = 0)=il-E Y. ^i 

\ j=[nsjj+l 



J2 {Zi - EZi) 

[ntj\ \ 



i{M^*j) = 0) - 1 + J2 ^i 



By equation (2.4) 
(2.5) 



k / [ntj\ \ k 

liml[[l-E Y. zA=l[{l-fiAAtl-si]). 



j=i 



i=[nsj\+l 



J=l 



By Lemma 2.1, 



(2.6) 



E 



fl-l \ / L"*d 

Vi=l / Vi=Lns;J+l 



as n — > oo. Now, Yl,iL\ns\+i^j=i+iZiZj is the number of pairs of points 

in i/„(A*) and so is equal to ^i^„(Ap(i/„(Ap — 1). This is greater than or 
equal to f„(A^) — 1 when f„(A*) > 1, so 

\ntj\ L'^ijJ L'^^jJ 

0</(i/„(A*) = 0)-l+ E ^^^ E E ^*^i 

i=\nsj\+l i=[nsjj+l j=*+l 
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and so, by Lemma 2.1, 



limsupE' 

n 

(2.7) 



\j=l ) \ j= Ins, 1+1 / 



<\A\ii-sif. 



Putting together equations (2.5), (2.6) and (2.7), we obtain 

k 
E[VI{MA*) = 0)] -EVUil- /.A, [ti - si]) 



lim sup 



<e. 






Taking increasingly fine partitions of A*, Hj^iCl ~ ^aJ^; — si]) = Y\j^i{l — 
i/'(Ap) converges to exp(— i/'(Ap) = EI{v{A*) = 0), wliich completes the 
proof. D 

Now, if A* is a finite union of finite rectangles in (R\(— c, c)) x M, then 
it can be written as A* = [J^^lAi x [si,ti), where Aj are finite unions of 
finite intervals in ]R\(— c, c) and si < ti < S2 < ■ ■ ■ ^ ^fc- Then, by Lemma 2.1, 
Ei^niA*) -^ Ya^iEv{A*) and, by Proposition 2.1, EI{vn{A*) = 0) ^ 
EI{i/{A*) = 0) as n^oo. Also, lima;_oolimsup„P(fn(M\(-rc,a;) x [0,t])) = 
0, so, by [9] Theorems 4.7 and 4.9, f„ converges weakly to u. 

Now, suppose that V is the indicator function of some event generated by 
Rn,H. Then, by Lemma 2.1, EVun{A*) -^ X;Li EVEu{A*) and, by Propo- 
sition 2.1, P(f„(A*) = 0) — > P(i/(A*) =0) as n ^> oo. So, asymptotically, 
this convergence takes place independently of V and we can conclude that 

(2.8) (L(n)-''n-(i-'^(i-^))5]S(/(X,)|.F.-d),z.J ^(/,i2.,^(t),i/) 

weakly, jointly where R^^h and i' are independent. 
By Karamata's lemma (see [10]), 

EifiX,)Ii\fiX,)\ < can)f ~ '^c^-^aln-^ 

Z — a 

as n ^ cx). If Yi = f{X,)I{\f{X,)\ < can) - Eif{Xi)I{\f{X,)\ < ca„)|.Fi„rf), 
then Yfc, Ifc+j^, Yk+2d, ■ ■ ■ are martingale differences for each 0< k < d. Then, 
using Doob's maximal inequality, 

/ [nt] x2 _ 

(2.9) limsup^ maxa-^VyJ <4dt'^- c^"". 

^ ' n \o<t<T " ^ 7 - 2- a 

Now, define ¥( = f[Xi)I{\f{Xi)\ > can)- Again, by Karamata's lemma, 

E\Yir'r.-^^^canf^<>--L,{can) 
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as n — > oo. Applying Theorem 5.1 of [8], we have 

/ LntJ \ 2 

Eia-'Y.E(^^\^^-d)-EYn 

< C(ni)2^Li(nt)2a;;2a/(i+e)^^(^^^)2/(i+9) ^ q^ 
which imphes tightness and finite-dimensional convergence, so 

[nt\ 

(2.10) a-' Y. iE{Yl\r,.d) - EYl) ^ 

i=l 

weakly. 

Now, noting that 

lim limsupP(i/„(IR\(-x,x) x [0,tl)) = 0, 

X — ^oo ^ 

it follows that 



(2.11) a-i Y. /(X.)/(a-i|/(^.)l > c) ^ / ydu 



weakly in the Skorohod topology. By the representation given in [11], Sec- 
tion 3.12, 

(2.12) f fyd{u-u')^R*{t) 
Jm\{-c,c) Jo 

in probability in the Skorohod Ji -topology as c — > oo. Applying Theorem 4.2 
of [1] to equations (2.9), (2.10), (2.11) and (2.12), we conclude that 

[ntj 

(2.13) a-' Y fi^i) - E{f{X,)\Ti_d) ^ R*{t) 

i=l 

weakly as n — > oo. Because of the asymptotic independence of i?K,i-K(i-H) 
and I', 

( \nt\ \nt\ \ 

V 1=1 j=i / 

where i?K,i-K(i-_f/)(i) and R*{i) are independent. Theorem 1.1 follows im- 
mediately. 
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The proof of Theorem 1.2 proceeds similarly. The measure Un can be 
constructed as before. To apply Theorem 4.2 of [1] and complete the proof, 
we need only show that 

/ Int] \ 2 

(2.14) limlimsup E i max a~^y^Yi-EYi] =0, 

^ ^ c^o n \o<t<T {r[ I 

where Yi = f{Xi)I{\f[Xi)\ < can), which follows from showing that the lower 
order chaos terms are insignificant in the limit. 
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